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ABSTRACT 

Recent observational work by Israelian et al. has shown that sun-like planet host stars in the temperature range 5700 K< Tgff < 
5850 K have lithium abundances that are significantly lower than those observed for "single" field stars. In this letter we use stellar 
evolutionary models to show that differences in stellar mass and age are not responsible for the observed correlation. This result, along 
with the finding of Israelian et al., strongly suggest that the observed lithium difference is likely linked to some process related to the 
formation and evolution of planetary systems. 
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1. Introduction 

The increasing number of discovered planets continues to 
give new insights i nto the planet formation processes (e.g. 
lUdrv & Santosl2007h . The study of planet host stars is providing 
crucial observational evidence for this. The first large uniform 
studies (Santos et al. 2001, 2004b) have shown that stars with gi- 
ant planets have chemical abundances that are distinctly different 
from those found in "single" stars, a resu lt that was later con- 
firmed by other authors (e.g. iFischer & V alenti 2005). This re- 
sult provided important constraints for the mod els of planet for- 
mation and evolu tion (e.g. Pollack et al.. 1996t iMordasini et al.l 
|2009HBossII2001 . 

Although most studies of chemical abundances in stars with 
planets have concentrated on the measurement of iron abun- 
dances as a metallicity proxy, several works have been made to 
study the abundances for a vari ety of refractory and volatile ele- 
ments in planet host stars (e.g. Gilli et all 2006t lEcuvillon et al.l 
l2006tlTakeda et al.l2007l:lNeves et al.l2009l) . Besides the general 
chemical enrichment found for stars hosting giant planets (inter- 
estingly not found for stars hosting very low mass planets, see: 
ISousa et al.ll2008h . none of these found compelling evidence for 
other chemical peculiarities. 

The possibility that stars with planets present different abun- 
dances of the light elements lithium and beryllium has also 
been debated ([King et al. 1997t Garcia Lo pez & Perez de Taoro 
198'; 'Delivannis et al."2000'; 'Cochran et al. 1997; Ryan 2000; 
,Gonzalez & Laws 2000; Gonzalez 2008; Israelian et al. 2004; 
Santos et al. 2004a; Takeda & Kawanomoto 2005; Taked a et all 
2007;,Chen & Zhaa2006:Luck & Heiter 2006 ). Different abun- 



dances could indicate that planets or planetary material were en- 
gulfed by the star during its lifetime (and in which quantity) (e.g. 
israelian et al. 2001), or suggest that the rotational history of the 
stars depends on the existe nce of planets or indirectly on the pro- 
cess for planet formation (lBouvierl|2008t ICastro et al.ll2009l) . 



Recently, llsraelian et al.l (l2009h presented a large uniform 
study of lithium abundances in a sample of stars from the 
HARPS planet search programme dMavor et al.ll2003h . In their 
work they reported a significant difference regarding the de- 
pletion of lithium for planet host stars when compared with 
stars with no detected planet in the ran ge 5700 K< Tf ff < 
5850 K, c onfirming former suspicions (e.g. llsraelian et alll2004l; 
iTakeda et al. 2007). According to Israelian et al., stars with plan- 
ets in the temperature range around the solar temperature (solar 
analogues) have significantly lower lithium abundances when 
compared with "single" stars (for which no planets were de- 
tected so far). The uniformity of the HARPS sample (composed 
mainly of old inactive stars) allowed Israelian et al. to exclude ef- 
fects like stellar rotation, stellar activity, or chemical abundances 
as the cause for the observed difference. 

In this paper we use stellar evolution models to explore the 
possibility that stellar age and mass could be responsible for the 
observed difference. In Sect. 2 we present the procedure used for 
the determination of precise and uniform masses and ages for 
our sample. In Sect. 3 we explore possible correlations between 
these two parameters and the lithium abundances. We conclude 
in Sect. 4. 

2. Masses and ages 

The stellar age and mass was determined by means of com- 
parison between stellar evolutionary models and observations. 
Individual values of the stellar luminosity, effective temper- 
ature, and metallic ity for the solar analogues in our sample 
('israelian et al.' '2009) were taken from the uniform study of 
(Sousa et al. 2008). These were used for a comparison wi th stel- 
lar models computed with the CESAM code version 3 (iMorell 
1997)j running in the Coimbra Observator y. The details on the 
physic s of these models can be found in iFernandes & SantosI 
(l2004l) with one exception: the equation of state FEE is used 
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(lEggleton et alJ [T973I) . an analytical equation of state suitable 
for FGK stars. With these input physics the solar model fits the 
observed luminosity and effe ctive temperature for the c ommon 
accepted solar age of 4.6 Gyr dPziembowski et al.lll999h . as bet- 
ter as 10 considering a mixing length parameter of convection 
a=1.65, the helium abundance 7=0.2675, and the metal abun- 
dance Z=0.0173. With these values, and assuming the primor- 
dial helium equal to 0.25 (e.g. lPeirnbe rt et al. 2009), the helium 
abundance to metallicity ratio is A F/A Z = 1.53. The mass and 
age resolution for the models used in this procedure are 0.01 Mq 
and a maximum time step of 200 Myr respectively. 

For each star with a metal abundance fixed to the observed 
metallicity value (according to Sousa et al. 2008), we computed 
several models for different ages and masses to fit the observed 
luminosity and effective temperature. We assumed for those 
models the solar mixing length parameter, a null overshooting 
parameter (aov) for stars with masses lower than 1.1 Mq, and 
Oov - 0.2 for higher mass stars (e.g. Claret.200 7). 

The stellar mass and age are estimated as follows: for each 
stellar evolutionary model we compute the functional 



log Lobs - log L„, 



where 
values 



O-L 

the subscripts obs 



log Teffobs - log Teff„„d 



Teff 

and mod refer to the observed 
(model) values respectively (e.g. 
Lastennet & Valls-Gabaudll2002 ). The cr values are taken from 
Sousa et al. I (120081) . The solution is obtained by minimizing the 



and theoretical 



above function. For each star the final mass and age is obtained 
when the prediction of the theoretical stellar model is inside the 
observational error bar, \.t.,x^{M, t) < 2. 

In Table [T] we present the derived values for mass and age 
in our sample. The typical (relative) errors bars of our estima- 
tions are O.OIM© in mass and 0.5 Gyr and 1.0 Gyr, respectively, 
for ages lower and higher than 5 Gyr. This level of precision is 
possible because of the very small (relative) errors in the stel- 
lar pa rameters derived for the solar analogues in ISousa et aP 
(I2008h . We note however that we are not taking into account 
variations in parameters like the mixing length and helium abun- 
dance. Changes in these pa rameters could lead to erro rs in the 
stellar masses and ages (e.g. lFemandes & Santosll2b04l) . 

In order to check our methodology, we applied the stellar 
evolutionary models from the Padova group, computed with the 
web interface dealing with stellar isochrones and their deriva- 
tives (http://stev.oapd.inaf.it/cgi-bin/cmd) to the stars of our 
sample. Taking into account that the Padova models are con- 
strained to metal abundance lower than Z < 0.03, we have ar- 
bitrarily chosen a subsample of 27 stars to cover the range of 
mass and age of the stars in the paper In Figs. [T] and |2] we 
show the comparison between masses and ages respectively, us- 
ing the CESAM and Padova models. In addition we also inc lude 
a co mparison to the values published in 'Sousa et al.l (l2008l) and 
IValenti & Fischer (2005) (with this latter we have 41 stars in 
common). A very good ag reement exists between a ll sets of val- 
ues, except for the ages of IValenti & Fisched (l2005h which seem 
to be constrained to values closer to the Sun age. The reason 
for this is not clear; this discussion is out of the scope of the 
present paper However, if the age span of the stars in our sam- 
ple were smaller (taking the values from Valenti & Fischer), it 
would strengthen the results discussed below. 

3. Lithium vs. age 

Table[T] shows the values of mass and ag e derived for the star s 
in this work and their lithium abundances (ilsraelian et al.ll2009l) . 




1.0 1.2 
Mass (CESAM) [Mj 



Fig. 1. Comparison between masses de termined in this wo rk us- 
ing CESAM with the ones presented in lSousa et al. I dMl) (bot- 
tom), Valenti & Fischer (2005) (midle) and determined using the 
Padova models (top). The filled lines represent the identity line. 
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Fig. 2. Same as FigH] but for ages. 



The sample studied here corresponds to t he majority of solar 
analogues studied in Ilsraelian et al.l (l2009l) . However, here we 
only considered those belonging to the HARPS sample (~60). 
For these stars w e have the best uniform stellar parameters 
(ISousa et al ] I2008h . assuring the total consistency of our re- 
sults. Four stars in the sample were also excluded since no re- 
liable masses and ages could be derived (HD 28701, HD 78558, 
HD 96423, HD 143114). 

The top panel in Fig. |3] shows the lithium abundance plot- 
ted against the effective temperature for the stars presented in 
Table [U This r esult is almost indistinguishable from Fig. 1 in 
Ilsraelian et all (l2009l) . There is a general tendency for lower 
temperature stars to present on average lower lithium abun- 
dances, an expected r esult (see e.g. Pinsonneault et al. 199^ 
ISestito & Randichl2 005). It is also evident that the stars that host 
planets are distributed in the lower region of the plot, while the 
stars with no evidence for planetary companions are spread over 
the entire diagram. 

In the middle plot of Fig.|3]we present the lithium abundance 
as a function of the derived stellar age. As expected, the bulk of 
the stars are old. The HARPS sample was chosen to avoid active 
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Table 1. New mass and age values together with published values for lithium abundance (< represents upper limit), [Fe/H] and TefF. 
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(younger) stars due to th e higher difficulty i n searching for plan- 
ets in these objects (e.g. lSantos et al.ir2000l) . Interestingly a few 
younger objects are seen in the left part of this panel (<2Gyr), 
a region where no lithium-depleted star is seen. This suggests 
a lithium abundance dependence with age. Above an age of 
~2Gyr, however, we observe a wider dispersion in lithium for 
all stars in the plot, and no correlation is seen. 

The middle plot of Fig. 2 also clearly shows that even though 
we have planet hosts stars in all the age interval, these are prefer- 
ably located at lower lithium abundances. No particular age cor- 
relation is present. Age does not seem to explain the observed 
difference between the lithium abundances in the two samples. 

Finally, we present in the bottom panel the lithium abun- 
dance as a function the derived stellar mass. As for the age, we 
can easily see that the dispersion of the lithium abundances has 
the same behaviour for the complete range of masses in the sam- 
ple. No evidence is found that stellar mass could be responsible 
for the observed lithium offset between the two samples. We note 
that the two planet hosts with the highest lithium abundances 
show no particular age or mass tendency. Their high abundance 
is likely explained because they are among the highest tempera- 
ture stars in our sample. 

We have also verified if the relatively large range in metal- 
licities for our stars together with the dispersion in ages and 
masses could not be hiding any trend. To do this, we selected 
only stars within the metallicity interval 0.1<[Fe/H]<0.3, where 
the largest overlap between stars with and without planets exists. 
From these we selected only those with masses between 1.05 
and 1.15 Mo (Fig. 3, upper panel, lithium as a function of stel- 
lar age), and with ages between 2Gyr and 6Gyr (Fig. 3, lower 



panel, lithium as a function of stellar mass). The plots of Fig|4] 
do not show any hint of a mass or age-lithium correlation that 
could explain the observed difference between planet host and 
"single" stars. 

4. Conclusions 

We derived uniform values for the stellar ages and masses in a 
sample of solar analogue stars, for which Israe lian et aP (l2009h 
have found a clear difference in the lithium abundances corre- 
lated to the existence of planets. Our analysis has shown that age 
and mass cannot explain the lithium abundance differences. This 
strongly suggests that the observed differences are not related to 
stellar intrinsic properties. This result confirms the uniformity of 
the studied sample. 

Our results are directly linked with Li observations in open 
clusters. A large high-quality database for Li in open clusters ac- 
quired in the last years has permitted us to draw a secure picture 
of the empirical behav i our and evolut i on of Li in solar analogue 
stars jRandichl l2008t iRandich et all 120091; ISestito & RandichI 
l2005h . Numerous observations show that stars with the same 
age, temperature, and metallicity can be affected by different 
amounts of Li depletion. A dispersion in Li of at least a factor of 
10 has b een observed in the solar-age, solar metallicity clusters 
M 67 (P asquiniet ani2008h and several ot her old, s olar metal- 
licity or metal rich clusters (iRand ich 2008: Randich e t al. 200Sj). 
This spread strongly suggests that Li depletion must be affected 
by an additional parameter besides mass, age, and chemical 
composition. This parameter could be an initial angular momen- 
tum of the star affected by the formation and evolution of planets. 
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Fig. 3. Lithium abundances as a function of effective tempera- 
ture (top), age (middle), and mass (lower panel). Filled circles 
denote stars with planets, while open triangles denote "single" 
stars. AiTows indicate upper limits. 
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Fig. 4. Lithium abundances as a function of stellar mass (top) 
and age (bottom) in restricted subsamples of our stars. 



Similar Li dispersion (or often cited as bi-modality) is also ob- 
served in solar type fie l d star s (Favataetal. 1996; Galee v et al.l 
l2004t IChen et al.ll200Tl [Lambert & Reddvll2004 ). We think that 
the same mechanism is responsible for the large Li dispersion 
in the field and open cluster solar type stars. Figure |3] clearly 
supports this suggestion. 

No clear explanation has been found for the lithium abun- 
dance difference observed between stars with and without plan- 
ets, though a few possibilities have been suggested. These in- 



clude the star-planet interaction (ICastro et al.l l2009l) . the infall 
of planets into the star (leading to higher mixing - Theado et 
al., priv. comm.), or a difference in the rotational history of 
the star due to star-disc interaction. Massive proto-planetary 
discs capable of forming planets will likely help to break stel- 
lar rotation, thus changing the depletion rate of lithiurn (e.g . 
ICharbonnel & Primas"2005; Coc hran et all 1997 1 'Bouviei jl200a 
[Pinsonneault 2010). Interestingly, the same mechanism could 
have been responsible for the low lithium abundance observed 
in our own Sun, itself a planet host stars, and can be used to ex- 
plain the dispersi on of lithium abundance s observed in clusters 
of different ages dSestito & Randichll2005h . 
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